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Abstract 



The transverse energy (Ex) in PbPb collisions at 2.76 TeV nucleon-nucleon center-of- 
mass energy i^/s^) has been measured over a broad range of pseudorapidity (//) and 
collision centrality using the CMS detector at the LHC. The transverse energy den- 
sity per unit pseudorapidity (dEj/ dtj) increases faster with collision energy than the 
charged particle multiplicity. This implies that the mean energy per particle is in- 
creasing with collision energy. At all pseudorapidities the transverse energy per par- 
ticipating nucleon increases with the centrality of the collision. The ratio of transverse 
energy per unit pseudorapidity in peripheral to central collisions varies significantly 
as the pseudorapidity increases from f/ = to |//| = 5.0. For the 5% most central 
collisions the energy density per unit volume is estimated to be about 14 GeV/ fm"' at 
a time of 1 fm/c after the collision. This is about 100 times larger than normal nuclear 
matter density and a factor of 2.6 times higher than the energy density reported at 
^/s^ = 200 GeV at RHIC. 
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The goal of relativistic heavy-ion colHsions is to study the behavior of quarks and gluons under 
extreme conditions of pressure, density, and temperature, such as those that existed shortly af- 
ter the Big Bang. Similar conditions can be reproduced in the laboratory by colliding heavy nu- 
clei at the highest possible energies. Experiments at the Relativistic Heavy Ion Collider (RHIC) 
have shown that at nucleon-nucleon center-of-mass energy of y/s^ = 200 GeV a strongly in- 
teracting medium is produced. This system behaves as an almost perfect quantum fluid [1 -4J. 
There are also indications from the RHIC experiments that at these energies the initial state of 
the colliding nuclei may be a color glass condensate lllEllll. The presence of such a state may 
affect the spatial distribution of partons within the nucleus, particularly at high pseudorapid- 
ity [T, 'S]. Measuring the distribution of transverse energy over a wide pseudorapidity range 
sheds light also on the longitudinal expansion of the system. In 2010, the Large Hadron Collider 
(LHC) accelerated heavy ions to energies 14 times higher than RHIC in order to produce matter 
at energy densities never achieved before. One of the basic measurements in this new regime 
is that of the energy distribution of all the produced particles, which is connected to the initial 
energy and entropy densities of the produced matter. At lower energies, the measured rapidity 
distributions of particles are generally well described by Gaussians. The widths of these distri- 
butions are consistent with the predictions of Landau hydrodynamics, i.e. dy = ^y]nJ, where 
7 is the Lorentz factor of the colliding beams ||9l- [T3| . The rapidity variable, y = tanh^^(cz/c), 
where Vz is the velocity of the particle along the beam direction z, provides a way to describe 
the longitudinal distribution of matter created in these collisions. Calorimeters measure only 
the energy deposited at various angles and therefore the data are presented in terms of the dis- 
tribution of energy in pseudorapidity, rj. Pseudorapidity is defined as tj = — ln[tan(0/2)], with 
9 the polar angle with respect to the z axis. When the momentum of a particle is larger than its 
mass, its ?/ ~ J/. 

The Compact Muon Solenoid (CMS) experiment is a general-purpose detector designed to 
study hadron collisions at theTeV scale IIT4ll . In particular, it has almost hermetic calorime- 
try that is sensitive to the distribution of energy over nearly the complete angular range. The 
transverse energy is defined by Ej = E, sin 6i, where E, is the energy seen by the calorimeter 
for the i^'^ particle, 0, is the polar angle of particle i, and the sum is over all particles emitted 
into a fixed solid angle in an event. The quantity dEj/ dtj is an approximately Lorentz invari- 
ant measure of the energy distribution. The transverse energy is studied as a function of the 
geometry of the collision, i.e. the centrality, of the heavy-ion interaction. Finally, comparisons 
are made with lower energy data and theoretical models. 

The central feature of the CMS apparatus is a superconducting solenoid, of 6 m internal diam- 
eter, providing a magnetic field of 3.8 T. Within the central field volume are the silicon pixel 
and strip trackers, lead-tungstate crystal electromagnetic calorimeter (ECAL) and the brass- 
scintilla tor hadron calorimeter (HCAL). These calorimeters are physically divided into the bar- 
rel and endcap regions covering together the region of < 3.0. The Hadronic Forward (HE) 
calorimeters cover \rj \ from 2.9 to 5.2. The HE calorimeters use quartz fibers embedded within 
a steel absorber. The CMS tracking system, located inside the calorimeter, consists of pixel and 
silicon-strip layers covering |?/| < 2.5. A set of scintillator tiles, the Beam Scintillator Counters 
(BSC), are mounted on the inner side of the HE calorimeters to trigger on heavy-ion collisions 
and reject beam-halo interactions. In addition, two Zero Degree Calorimeters (ZDC) are used 
for systematic checks. Eor more details on CMS see |[T4|. 

In 2010, CMS recorded PbPb collision data corresponding to an integrated luminosity of 7.36 }ih' 
of which 0.31 jih^^ were included in this analysis. This luminosity selection provided a data 
sample with negligible statistical uncertainties. Minimum bias inelastic PbPb collisions were 
selected by requiring that either the HE or BSC counters detected a signal on both sides of the 
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interaction point. In the analysis at least two reconstructed tracks were required to form a ver- 
tex within ±25 cm of the nominal interaction point along the beam line and within a radius 
of 2 cm measured perpendicular to the beam relative to the average vertex position. Large- 
multiplicity beam-background events were removed by requiring the compatibility of the ob- 
served pixel-cluster lengths with the hypothesis of a PbPb interaction at the estimated vertex. 
Finally, events containing beam-halo muons were eliminated by a timing requirement on the 
BSC counters on opposite sides of the interaction point. The total event selection efficiency of 
the minimum bias trigger for hadronic PbPb interactions was found to be (97 ± 3)% [15J. 

Events were sorted into different centrality classes. The centrality of heavy-ion interactions is 
related to the number of participating nucleons and hence to the energy released in the colli- 
sions. In CMS, the centrality is defined as percentiles of the energy deposited in the HP. The 
most central /peripheral event class, i.e. (0-2.5)%/ (70-80)% in this analysis, has a large/small 
number of participants and a large /small energy deposit in HP. In order to estimate the mean 
number of participating nucleons ((Npart)) and its systematic uncertainty for each centrality 
class, a Glauber model of the nuclear collision was used IT6HT81I . 

The data were corrected for detector acceptance and inefficiencies using correction factors 
C(|?/|) estimated from the hydjet 1.8 |19| Monte Carlo (MC) event generator coupled to a 
Geant4 [20J CMS detector simulation. These correction factors were calculated as the ratio of 
MC predictions at the particle level and the detector level for each centrality class. The correc- 
tion factor C(|?/|) ~ 1.6 for \r]\ < 2, falls to ~ 1.1 by |//| = 4 and then rises to 2 at |7/| ~ 5. The 
non-linearity of the calorimeter response and the effect of the magnetic field cause the C ( | 1) 
to depend upon the pj spectra, the ratio of charged and neutral particles, and the mixture of 
mesons and baryons. The value of C(|7/|) increases if the assumed spectra shifts to lower pj 
or if the ratio of charged to neutral particles is larger. To estimate the systematic uncertainties 
in C(|//|), two tunes of HYDJET (1.6 and 1.8) were used, hydjet 1.8, (pT)ch = 0.66GeV/c, was 
tuned to LHC spectra and particle yields as measured by the ALICE collaboration [21J and 
successfully tested against a wide range of RHIC data, hydjet 1.6, (pT)ch = 0.57GeV/c, was 
only tuned to RHIC data llT9ll . At central pseudorapidity the fraction of Ej carried by charged 
pions, kaons, protons and antiprotons is 0.60 for HYDJET 1.6 and 0.62 for HYDJET 1.8. The re- 
sults were cross checked using data taken with no magnetic field and in addition, for B = 3.8 T, 
data from tracks with pj > 900 MeV/c were combined with energy clusters in the calorimeters 
to identify different types of particles and measure their energy. Since muons and neutrinos 
carry a negligible fraction of the total transverse energy and deposit almost no signal in the 
calorimeters they are not considered in this analysis and no correction factors are applied to 
account for them. The corrected transverse energy for N analyzed events is obtained as 

^(l^l) = c(M)- ^^(,,^^) (1) 

where the sum over / covers all calorimeters cells located within the rj range Arj and Ej^j is the 
transverse energy measured in a particular cell /. Note that for this summation no threshold 
is applied to the individual calorimeters cells. Several sources of systematic uncertainties were 
studied and their effects are summarized in Table [l] and described below. Energy Scale: all 
of the calorimeters were initially calibrated with test beam data and radioactive sources. Por 
the barrel and inner endcap calorimeters these calibrations were refined using isolated charged 
hadrons whose momentum was reconstructed in the tracker. Por the HP calorimeter the en- 
ergy scale was cross-checked by reconstructing Z — )■ e+e^ in pp collisions where either the 
positron or the electron was recorded in the ECAL and tracker. Symmetry about rj: for PbPb 
collisions the corrected dEj/ drj should be symmetric about rj = 0. The values of dEj/ drj for 
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positive and negative ?/ were found to differ by at most 0.5%. This close agreement implies 
that one can use the average of the two results as a best estimate of dEj/drj. Vertex distri- 
bution: the z distribution of the vertices is Gaussian with az ~ 6.1 cm. To test the sensitivity 
of Ej to the position of the interaction vertex along the beam line, z, the dataset was divided 
into two samples with \z\ < 10cm and 10cm < |z| < 25cm, respectively. The Ej distribu- 
tions of the two samples differ by less than 2%. Auto-correlations: since HF is used both to 
calculate centrality and to measure Ej for each centrality class, there is an auto-correlation in 
the measurement. This effect was estimated to be less than 1.5% by using a combination of the 
ZDC and pixel detectors to measure centrality. Calorimeter noise: the Geant4 simulation of 
the calorimeters included electronic noise. This noise was measured by studying a sample of 
events where the trigger only required the presence of clockwise and anticlockwise bunches of 
lead ions simultaneously in CMS. The simulation of the noise was checked by comparing the 
data to the simulated signal from a GEANT4 simulation of the most peripheral events in the 
data set. Any discrepancy in the simulation of the noise corresponds to a corrected average 
Ej per event of less than 5.8 GeV for \f]\ < 2.65 and 1.2 GeV for 2.65 < \t]\ < 5.2. This is only 
significant compared to the signal for (Npart) < 30. HF MC description: takes into account 
different ways of describing the dead areas of the HF detector. Centrality determination: the 
systematic uncertainty related to the centrality determination is only applied to the results that 
are normalized by (Npart)- 

Table 1: Systematic uncertainties for two regions in % for the most central (0-2.5)% 
((Npart) = 394) and most peripheral (70-80)% ((Npart) = 16) collisions. 
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Calo. noise 


14-18 


0.3 
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0.1-0.2 


HF MC 






1.5-9 


1.5-9 


Centrality 


6.7 


0.5 


6.7 


0.5 


Total 


14-22 


3.5-5.9 


11-17 


10-14 



Figure [T] shows the | rj \ dependence of the transverse energy density for four selected ranges 
of centrality. For the most central collisions ((Npart) = 394) dEj/dr] reaches 2.1 TeV at ?/ = 0. 
This is much larger than the value of 0.61 TeV measured at ^/s^ = 200 GeV |22| . At lower 
center-of-mass energies the pion multiplicity distributions are reasonably well described by 
Gaussians in rapidity with widths that are consistent with Landau-Carruthers hydrodynam- 
ics ||23l |24l. Since the mean pj of all particle species depends only weakly on rapidity, this 
implies that dEj/ dy is roughly Gaussian in rapidity at ^/s^ = 200 GeV. Recently, Wong has 
improved the formulation of Landau hydrodynamics [25] . This new formulation gives a better 
description of the 200 GeV RHIC data. At y/s;^ = 2.76 TeV and for {rjl < 5.2 the dEj/drj is 
consistent with a Gaussian (black solid line) with cr^ = 3.4 ± 0.1 for the most central collisions. 
The Gaussian and Landau curves in Fig. [l] are normalized to the CMS data at = 0. Both the 
Landau-Carruthers (blue dashed), and Landau-Wong (green dotted) formulations have distri- 
butions that are narrower than the data. Therefore the longitudinal expansion of the system is 
stronger than that predicted from either model. HYDJET 1.8, shown by the purple dashed line, 
has been tuned to LHC data in the small \f]\ region. It gives a good description of dEj/ drj at 
small I }] I but overestimates the data at large | r] \ for central collisions. The AMPT (A Multi Phase 
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Transport) model |l26ll27l (orange dashed line) overestimates dEj/ drj for central collisions but 
is in rough agreement with the shape of dEj/drj. Integrating [dEj/df])/ ((Npart)/2) over rj 
between —5.2 and 5.2 gives a total measured £t per participant pair of 80 ± 4GeV for the most 
central events. This serves as a lower limit for the total transverse energy per nucleon pair. 
Extrapolating to the full phase space gives a total transverse energy per pair of participating 
nucleons of 91 ± 5 GeV for the most central events. It is clear from Fig.jljthat the magnitude of 
dEj/df] increases rapidly with the number of nucleons participating in the collision. 
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Figure 1: Transverse energy density versus \rj\ distribution for a range of centralities of (0- 
2.5)%, (20-30)%, (50-60)% and (70-80)%. The boxes show the total systematic uncertainties. 
The statistical uncertainties are negligible. Also shown are a Gaussian fit and the predictions 
of various models (see text). The AMPT events are for perfectly central collisions. 

Figure |2] shows the evolution of (d£T/dJ/)/((Npart)/2) with (Npart) for several \r]\ regions. At 
all \rj\ values {dEj/ drj)/ ( (Npart) /2) increases with (Npart). This figure shows that the (Npart) 
dependence of transverse energy density changes as a function of pseudorapidity. This effect 
can be quantified by comparing peripheral (60-70)% ((Npart) = 30) to central (0-2.5)% colli- 
sions ((Npart) = 394) at various pseudorapidities. The ratio of peripheral to central (dEj/drj)/ ((Npart)/2) 
changes from 54 ± 2% at = to 68 ± 2% at\r]\= 5.0. The PHENIX collaboration at RHIC has 
studied transverse energy density in AuAu collisions for \ f]\ < 0.35 over a wide range of cen- 
tralities and for ^/s^ from 19.6 GeV to 200 GeV |l22l. At ^/s^^ = 19.6 GeV (dEj/d?/) / ((Npart) /2) 
at f/ = increases by a factor of 1.25 ± 0.17 as (Npart) increases from 63.8 to 336. At i/s^ = 
2.76 TeV this factor is found to be 1.47 ± 0.13 for a similar range of (Npart) • At ^/s^^ = 2.76 TeV, 
the HYDJET 1.8 code gives a good description of the centrality dependence of dEj/ drj ati] = 0. 

Figure |3] shows the energy dependence of {dEj/ drj)/ ((Npart)/2) for central collisions atrj = 0. 
For the top 5% most central events (dET/d)/)/( (Npart) /2) reaches 10.5 ± 0.5 GeV at a/s^ = 
2.76 TeV. The Ej rises more quickly with the center-of-mass energy than the logarithmic de- 
pendence used to describe data up to ^/s^ = 200 GeV [22J. For energies between 8.7 GeV 
and 2.76 TeV, dEj/drj at // = can be reproduced by a power-law dependence of the type 
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Figure 2: Transverse energy density normalized by ((Npart)/2) versus (Npart) for PbPb colli- 
sions at = 2.76 TeV at several values of |?/|. The bands show the total systematic uncer- 



tainties. The statistical uncertainties are negligible. For rj 
results from the HYDJET 1.8 model are also shown. 



lower-energy PHENIX data and 



0.2. A similar effect has been seen in the measurement of the ^s^^ evolution of 
the charged particle multiplicity IIT8l l28| . The {dEj/dtj)/ ( (Npart) /2) increases by a factor of 
3.07 ± 0.24 from ^s;;;;; = 200 GeV to 2.76 TeV. This is to be compared to a factor of 2.17 ± 0.15 
for the pseudorapidity density (dNch/d//)/ ((Npart) /2) ||l8l|2Tl|22l. For the 5% most central 
collisions CMS has measured dEi/d?/ = 2007 ± 100 GeV and dNch/df/ = 1612 ± 55 HI. Divid- 
ing the measured transverse energy by the observed charged particle multiplicity for the same 
centrality gives a transverse energy per charged particle of 1.25 ± 0.08 GeV at ^/s^ = 2.76 TeV. 
This compares to 0.88 ± 0.07 GeV at = 200 GeV [22J. 

The sum of the transverse energies of all particles produced in the event depends upon both 
the entropy and temperature of the system. Using geometrical considerations, Bjorken II37I 
suggested that the energy density per unit volume in nuclear collisions could be estimated 
from the energy density per unit rapidity. A commonly used estimate of energy density is 
given by l|22[ 

where A is the overlap area of the two nuclei and Tq is the formation time of the produced sys- 
tem. The Jacobian J{y,rj) depends on the momentum distributions of the produced particles. 
In the limit that the rest mass of the particles are much smaller than their momenta J{y,tj) = 1. 
The average Jacobian was calculated using hydjet 1.8 for \ f]\ < 0.35. For central collisions at 
^s;^ = 2.76 TeV, J{y,ri) = 1.09. This compares to 1.25 at = 200 GeV [22]. For the top 

5% most central collisions this formula gives e = 14 GeV/ fm at a time tq = 1 fm/ c and for a 
transverse surface of A = tt x (7fm)^ [|22J . This is a factor of 2.6 times larger than the energy 
density calculated at = 200 GeV Il22l . 

In summary, for the most central PbPb collisions at a/s^ = 2.76 TeV, the maximum of the 
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Figure 3: Transverse energy density normalized by ((Npart)/2) for central collisions at = 
versus a/s^. The vertical bar on the CMS point is the full systematic uncertainty. The statistical 
uncertainty is negligible. The solid line is a parametrization of lower-energy data compiled by 
the PHENIX Collaboration, see p2ll29ti36| . The dashed line corresponds to a power law fit s'^^ 
for ^s;;;^ > 8.7 GeV. 

transverse energy distribution has been found to be 2.1 TeV at // = 0. Even at a very forward 
pseudorapidity of \rj\ = 5.0, dEj/ drj and hence the energy density of the produced system at 
the LHC is larger than that measured for ?] = at RHIC. At 2.76 TeV, the shape of dEj/ drj is 
consistent with a Gaussian function of width a,^ = 3.4 ± 0.1 for central collisions. This distri- 
bution is wider than the prediction of Landau hydrodynamics but narrower than that given by 
the HYDJET 1.8 simulation. The (d£T/d?/)/((Npart)/2) increases with (Npart) at all pseudora- 
pidities. The ratio of transverse energy in peripheral compared to central collisions increases 
by a factor of 1.26 ± 0.06 from t] = to \rj\ = 5. The transverse energy density at ?/ = grows 
more rapidly with the center-of-mass energy than the logarithmic scaling with ^Js^ that de- 
scribes lower-energy data. It also grows faster with energy than the multiplicity, implying a 
significant increase of the mean transverse energy per particle compared to lower energy data. 
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